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Abstract. The model dependence of the development of extensive air showers 
generated by high-energy cosmic-ray particles in the atmosphere is studied. The 
increase of proton-proton and proton-air inelastic cross-sections and values for 
the elasticity are varied in the hadronic interaction model QGSJET. Using the 
CORSIKA simulation program, the impact of these changes is investigated on air 
shower observables like the average depth of the shower maximum Xmax and the 
number of muons and electrons at ground level. Calculating the mean logarithmic 
mass from experimental Xmax values, it is found that a moderate logarithmic 
increase of the proton-proton inelastic cross-section from a^"^ = 51 mb at 
Eo = 10^ GeV to o-^^^' = 64 mb at Eq = 10* GeV and an elasticity additionally 
increased by 10% to 15%, describes the data best. Using these parameters, the 
mean logarithmic mass (In A) derived from Xmax measurements is compatible 
with the extrapolations of the results of direct measurements to high energies 
using the poly-gonato model. 
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1. Introduction 

Cosmic rays above energies of 10^"* eV are investigated by measurements of extensive 
air showers in the atmosphere. These cascades of secondary particles are generated 
by interactions of high-energy particles with air nuclei. While the shower develops 
in the atmosphere the number of particles increases up to a depth where the average 
energy of the secondaries equals a critical energy. Beyond this point the number of 
particles decreases approximately exponentially. The depth where the cascade reaches 
the largest number of charged particles is usually referred to as the depth of the shower 
maximum Xmax- 

Two basic methods are used to examine air showers. In the first, the debris 
of the particle cascade on the Earth's surface is observed in large detector arrays, 
measuring electrons, muons, and hadrons, viz. particle distributions at ground level. 
The second, calorimetric method investigates the longitudinal development of the 
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shower by registration of Cerenkov or fluorescence light, generated by the shower 
particles in the atmosphere, and eventually X^ax is derived. 

Recently, the author has compiled the results of many air shower experiments with 
the intention to estimate the best primary all-particle spectrum (Horandel 2003). For 
this purpose, the energy spectra of individual elements as directly measured at the top 
of the atmosphere below 1 PeV have been extrapolated to high energy using power 
laws and assuming individual rigidity dependent cut-offs. The all-particle spectrum 
of these direct measurements was compared with air shower observations at higher 
energies. For the indirect measurements it turned out that the individual energy scales 
had to be renormalized to match the all-particle spectrum of direct measurements at 
1 PeV. Only small energy shifts were necessary, all within the errors quoted by the 
experiments. But most shifts had a negative sign, indicating an overestimation of the 
primary energy. A likely explanation for this effect will be discussed later in chapter 6. 

As a result all experiments yielded consistent all-particle energy spectra. The 
extrapolations of the direct measurements have been fitted to the average all-particle 
flux of the indirect measurements to determine the parameters of the individual 
cut-offs mentioned above. It has been shown, that the experimental data can be 
parametrized consistently within a phenomenological model, the poly-gonato model. 
The mean logarithmic mass (In A) obtained from the model is compatible with results 
from experiments measuring particle distributions at ground level. But the mass 
composition disagrees with results from experiments investigating the longitudinal 
shower development in the atmosphere. 

A possible cause for the discrepancy might be, that the measured cascades 
penetrate deeper into the atmosphere than predicted by the simulations, viz. the 
codes predict a too small value of Xmax- A similar conclusion was drawn by Erlykin 
and Wolfendale (2002). Investigating the muon production height with the KASCADE 
experiment, Biittner et al (2001) also find an indication for a deeper penetration of 
cascades into the atmosphere as compared with predictions of the air shower simulation 
program CORSIKA (Heck et al 1998) with the high-energy interaction model QGSJET 
(Kalmykov et al 1997). 

Reasons why the codes predict too small Xma^-values may be numerous. 
Objective of the present article is to study the influence of model parameters like the 
inelastic cross-sections or the elasticity on the average value of X^ax and, consequently, 
on the mean logarithmic mass derived from it. Goal of the investigations is to perhaps 
resolve the discrepancies in (In A) outlined above. 

Some aspects of the impact of the complex system of parameters in simulation 
codes on Xmax are briefly summarized in chapter 2. Variations of the interaction 
model QGSJET are described in chapter 3. The resulting changes of Xmax and the 
consequences on (hi A) are presented in chapters 4 and 5. For completeness, the 
influence of the changes on the number of electrons and muons at ground level is 
briefly sketched in chapter 6. 

2. Model peirameters and the average depth of the shower maximum 

The longitudinal development of the nuclear cascade depends among others essentially 
on following physics quantities: The inelastic cross-sections (Jmei of primary and 
secondary particles with air nuclei, the average number of particles produced in an 
interaction — the multiplicity /U, and the average fraction of energy transferred into 
secondary particles — the inelasticity K. An increase of the inelastic cross-sections 



Total inelastic cross- sections and X^ax in extensive air showers 



3 



results in an earlier development of the cascade. A similar effect has an increased 
inelasticity. The particles lose more energy, hence the shower reaches its maximum 
earlier in the atmosphere. A reduction of the multiplicity produces less particles in 
the first interactions. Due to energy conservation they will be more energetic and the 
shower develops more slowly, i.e. the maximum is deeper in the atmosphere. Pajares 
et al (2000) derived a parametrization for the relative change of Xmax as function of 
the changes of /x and K 

AXma^ _ lAfl 1 AK 
Xmax ^ 2 f, 10 K ■ ^ ' 

A similar relation for the dependence of Xmax on inelastic cross-sections will be given 
below in chapter 4. 

A variety of models has been used to study the development of extensive air 
showers and the effects of different parameters in the simulations on Xmax have been 
elaborated in the literature. Some of the findings are summarized in the following. 

The average depths of the shower maximum obtained for different interaction 
models implemented in CORSIKA — i.e. DPMJET 2.5 (Ranft 1995 and 1999), 
NEXUS 2 (Drcschcr et al 2001), QGSJET 01, and SIBYLL 2.1 (Fletcher et al 1994, 
Engel et al 1999) — have been compared by Knapp et al (2003). At Jcnee energies a 
spread of the models of about 50 g/cm^ has been found for proton induced showers 
and 20 g/cni^ for primary iron nuclei. 

Similar results have been obtained by Fowler et al (2001). The systematic 
differences in Xmax between different models in CORSIKA — HDPM (Capdevielle 
et al 1992), QGSJET, SIBYLL, VENUS (Werner 1993) — at 1 PeV are specified as 
45 g/cm^ for primary protons and 25 g/ cm^ for iron nuclei. 

The longitudinal shower development for the hadronic interaction models SIBYLL 
and QGSJET embedded in the framework of the MOCCA (HiUas 1997) and CORSIKA 
programs has been explored by Pryke (2001). At 1 PeV the differences in Xmax 
between CORSIKA/QGSJET and MOCCA with its internal hadronic interaction 
model amount to 45 g/cm^ for primary protons and to 30 g/cm^ for iron induced 
showers. The differences between the interaction models are related to distinct 
inelasticities: QGSJET produces more inelastic events, which lead to less deeply 
penetrating showers. Using the MOCCA frame in connection with its internal and 
the SIBYLL interaction model, the values differ only by about 5 to 10 g/cm^. This 
is compatible with estimates of Dickinson et al (1999) who obtained an increase of 
about 10 g/cm2 from MOCCA/MOCCA to MOCCA/SIBYLL. 

Wibig (2001) discussed the influence of various hadronic interaction models — 
FRITIOF (Anderson et al 1991), CMC (Wibig 1997), HDPM, QGSJET, SIBYLL, 
and VENUS — on the predictions for Xmax using the shower simulation program 
CORSIKA. For a primary energy of 1 PeV he found differences in the order of 40 g/ cm^ 
for primary protons and 18 g/cm^ for iron nuclei. 

A fast one-dimensional hybrid method has been used to simulate air showers 
by Alvarez-Muniz et al (2002a). The models QGSJET 98 and two versions of 
SIBYLL (1.7 and 2.1) have been applied to describe the hadronic interactions. For 
SIBYLL 2.1 the inelastic proton-air cross-section rises faster with energy as compared 
with QGSJET, whereas the inelasticities are almost equal in both models. On the 
other hand, the multiplicity of charged secondary particles produced in proton-air 
collisions grows rapidly with energy for QGSJET. Compared with QGSJET, for 
SIBYLL 2.1 the mean values of Xmax of proton induced showers are about 6 g/cm^ 
larger at energies of 1 PeV and 10 g/cm^ at 1 EeV. 
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Gaisser et al (1993) investigated the statistical model (Landau 1969, Fowler 
1987), the Kopchovich-Nikolaev-Potashnikova (KNP) model (Kopcliovich 1989), and 
the mini-jet model (Gaisser and Halzen 1987). The characteristic differences between 
them are the energy dependences of the inelasticity. In the statistical model the 
inelasticity decreases with energy, almost constant values arc obtained for the mini-jet 
model, and in the KNP model the inelasticity grows with energy. At 1.25 EeV the 
depth of the maximum of proton initiated showers varies by about 45 g/cw?. Due to 
the larger inelasticity the showers develop faster in the KNP model as compared with 
the statistical model. The latter represents the other extreme, generating relatively 
long showers. 

The impact of internal parameters in the Quark Gluon String model (QGSJET) 
on the average depth of the shower maximum was examined by Kalmykov et al 
(1995). Different implementations for the contributions of semihard processes lead 
to differences in Xmax for primary protons of about 18 g/cw? and 40 g/crc? at 1 PeV 
and 100 PeV, respectively. The role of the production of charm particles for the air 
shower development was found to be negligible, the estimated differences are in the 
order of 3 g/cm^. 

The effect of the multiplicity of secondary particles on X^ax was scrutinized by 
Anchordoqui et al (1999) using the AIRES program (Sciutto 1998) with the models 
SIBYLL and QGSJET. For both models the AIRES cross-sections have been utiHzed 
in the calculations. As a consequence of the lower inelasticity in SIBYLL, the model 
produces fewer secondaries than QGSJET. Hence, there is a delay in the shower 
development for SIBYLL generated cascades. Primary protons with energies of 1 PeV 
penetrate deeper into the atmosphere by about 45 g/cm^ and 100 g/cw? at 1 EeV. 

Capdevielle and Attallah (1995) pointed out how uncertainties in the description 
of parton distribution functions of hadrons influence the longitudinal development of 
air showers. The effects on inelastic cross-sections, multiplicity and inelasticity were 
derived. The maximal uncertainty in X^ax was estimated to be at 1 EeV in the order 
of 95 g/cm^ for proton induced showers and 40 g/cm^ for iron primaries. 

A second step cascading mechanism was discussed in the geometrical multichain 
model by Wibig (1999). Second step cascading is defined as the interaction of a 
wounded nuclcon of one nucleus with another nuclcon from the same nuclciis before 
hadronization occurs. These internuclear cascades result in an increase of the hadron 
air inelasticity for heavy nuclei. Results of CORSIKA simulations with this interaction 
model indicate that this mechanism results in an earlier shower development for iron 
nuclei of 10 g/cw? to 30 g/cm^ at energies of 10 PeV to 10 EeV. 

New effects in high-energy interactions like percolation, quark gluon plasma, or 
string fusion were anticipated by Pajares et al (2000). The authors concluded that 
these effects dump the multiplicity and increase the inelasticity. This leads to larger 
values oiXmax- 

Erlykin and Wolfendale (2002) linked the observed discrepancies in (In A) to two 
effects in nucleus- nucleus collisions not been taken into account before: A few percent 
energy transfer into the electromagnetic component due to electron-positron pair 
production or electromagnetic radiation of quark gluon plasma and a small slow-down 
of the cascading process in its initial stages associated with the extended lifetime of 
excited nuclear fragments. The latter displaces the shower deeper into the atmosphere. 

To summarize, several hypothetical effects have been described which could 
change the longitudinal shower development and their influence is not easy to 
segregate from each other. The maximum model ambiguities, taking into account all 
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considerations as discussed, amount in the Jaiee region to about AXma^ ~ 45 g/cm^ 
for primary protons and 25 g/cm^ for iron induced showers. These imcertainties 
increase with energy to /^.X^ax ~ 100 g/cvc? for proton induced showers at 1 EeV. 
In principle, the quoted ambiguities could be reduced (by some 30%) as some 
model approaches can be rejected by theoretical arguments and/or comparison with 
accelerator or cosmic-ray data. In addition to the model ambiguities discussed so far, 
there are experimental uncertainties, which will be discussed in chapter 4. 

3. Inelastic hadronic cross-sections and other model parameters 

Among the various parameters which control the shower development, the inelastic 
cross-sections and the elasticity will be scrutinized more closely. The actual situation 
can be recapitulated as follows. 

3.1. Inelastic cross- sections 

Inelastic cross-sections for proton-proton interactions are of great interest for particle 
physics. They have been measured in detail in collider experiments, but studied 
also at highest energies using cosmic rays. At energies exceeding ^/s = 1 TeV 
the knowledge about the increase of the proton-proton cross-section as function of 
energy is limited by experimental errors of both, collider and cosmic-ray experiments. 
The insufficient knowledge of the high-energy cross-sections restricts the reliability of 
simulation programs to calculate the interactions of high-energy cosmic rays. 

A recent review of pp cross-sections from collider experiments has been given by 
Hagiwara et al (2002). At present, the highest energies at colliders are available at 
the Tevatron = 1.8 TcV). There, values for the total pp cross-section disagree 
by about 10% between different experiments, i.e. 71.71 ± 2.02 mb (E-811, Avila et al 
1999), 72.8 ± 3.1 mb (E-710, Amos et al 1992), and 80.03 ± 2.24 mb (CDF, Abe et al 
1994a). 

In cosmic-ray experiments the attenuation of the proton flux entering the 
atmosphere is used to derive the inelastic proton-air cross-section. Then, the latter is 
utilized to derive the proton-proton cross-section using the Glauber theory (Glauber 
and Matthiae 1970). In air shower experiments three basic methods are applied to 
obtain cross-sections: 

1) The ratio of the proton flux at the top of the atmosphere $o to the flux of surviving 
protons $g measured below the atmosphere at a depth x. From these quantities 
the interaction lengths of protons in air Xp^air{E) is determined using the relation 
^g{E,x)/^o{E) = exp[—x/ \p-air{E)]. This method needs high fluxes of penetrating 
primary protons and, therefore, can be applied at relatively low energies only, see 
Yodh et al (1983). 

The second and third method take advantage of the fact that cosmic-ray protons 
interact in the atmosphere at rates which decrease exponentially with increasing depth. 

2) For a fixed primary energy and zenith angle the distribution of the average depth 
of the shower maximum has an exponential tail with a slope given by the attenuation 
length A, see e.g. (Baltrusaitis et al 1984) or (Gaisser et al 1993). 

3) The last method to derive A utilizes the zenith angle distributions of the shower 
intensity for fixed primary energies, as described by Hara et al (1983). 
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Figure 1. Inelastic proton-air cross-sections 
versus primary energy. For references of the 
experimental values see text. The lines repre- 
sent calculated cross-sections for three different 
versions of the model QGSJET 01, model 1 

(original QGSJET) ( ), model 2 ( ), 

and model 3 ( ), see text. 



The attenuation length is connected with the interaction length A by the relation 
A{E) = k{E)X,_^UE) = HE) ^t:;^ . (2) 

k{E) is a model dependent proportional factor, which among others depends on the 

mean inelasticity of the interactions K{E). o-^iiriE) is the total inelastic proton-air 
cross-section, {Aair) = 14.5 the effective atomic weight of air with the proton mass 

The conversion of the attenuation lengths to cross-sections depends essentially 
on theoretical assumptions on k{E), for a discussion and theoretical motivation see 
e.g. (Kopeliovich et al 1989), (Engel et al 1998), or (Block et al 2000). For example, 
using a different value for k in equation (2), total inelastic cross-sections obtained by 
Baltrusaitis et al (1984) and Honda et al (1993) around cm. energies ^/s ~ 10 TeV 
have been reduced by 10% to 15% by Block et al. The original value of the Fly's 
Eye experiment cr^^^^ = 530 ± 66 mb at ^/s = 30 TeV has been rescaled to 
cji'^e].^ = 460 ± 40 mb. 

The subsequent conversion of proton-air to proton-proton cross-sections depends 
on theoretical uncertainties as well. The details required to apply Glauber theory 
in this context have been discussed by Gaisser et al (1987) as well as by Wibig and 
Sobczynska (1998). The latter apply a geometrical scaling hypothesis, use an exact 
Glauber formalism and conclude, that the proton-proton cross-sections reported by 
Baltrusaitis et al and Honda et al are overestimated by about 10%. 

Besides the theoretical uncertainties, all methods require that the primary cosmic- 
ray flux contains a sufficient fraction of protons, e.g. Hara et al (1983) derive cross- 
sections under the assumption of a contribution of at least 10% protons. Thus, 
additional systematic errors have to be taken into account concerning the unknown 
mass composition at high energies. 

In early cosmic-ray work an increase of the inelastic cross-section as function of 
energy has been found by Yodh et al (1972, 1973). In the energy range up to 100 TeV 
the measured flux of unaccompanied hadrons at ground level has been used to calculate 
lower bounds for the proton-air cross-sections by Yodh et al (1983) combining several 
experimental results and also by Mieike et al (1994) using the KASCADE prototype 
calorimeter. Particle numbers measured with the EAS-TOP experiment have been 
used to deduce a value for the inelastic proton-air cross-section at knee energies 
(Aglietta et al 1997). At the highest energies <tI^I^^^ has been derived from the shower 
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Figure 2. Relative fraction of protons in 
cosmic rays versus primary energy according 

to the poly-gonato model ( ). The band 

illustrates the expected error. The line above 3- 
10^ GeV is an upper limit assuming the ad-hoc 

component being protons only ( ). The 

points are obtained by fits to measured Xmax 
distributions of the Fly's Eye experiment, see 
section 4.2. They are connected by a line to 
guide the eye ( ). 

attenuation length A by the Fly's Eye experhiient (Baltrusaitis et al 1984) using the 
exponential tail of Xmax distributions and by the AGASA experiment (Honda et al 
1993) measuring particle numbers at ground level. Dyakonov et al (1990) derived lower 
limits for cross-sections from investigations of the Yakutsk data, applying two methods, 
fits to the exponential tail of Xmax and observing the zenith angle distribution of the 
shower intensity for fixed energies. A second analysis of the Yakutsk data (Knurenko 
et al 1999) takes into account the tail of the Xmax distributions at fixed energies. The 
results obtained are shown in figure 1 as function of laboratory energy Eq ■ 

The experimental values shown in figure 1 have been derived assuming a 
sufficiently large contribution of protons in primary cosmic rays. The validity of 
this assumption shall be examined next. For this reason, the fraction of protons as 
function of energy is shown in figure 2. It has been obtained from the analysis of direct 
measurements and air shower data applying the poly-gonato model, as outlined in the 
introduction. The solid line represents the ratio $p/$a/; for the cosmic rays which shall 
be denoted, with precaution, as the "galactic component". The error band indicates 
the uncertainties expected from the model. Since the average all-particle spectrum is 
reasonably well known, the error indicated depends essentially on the parametrization 
of the proton fiux and its errors. According to the poly-gonato model the galactic 
component is characterized by cut-offs of the flux of individual elements at energies 
proportional to their nuclear charge {Ez = Ep- Z), starting with protons at an energy 
Ep — 4.5 PeV. Above 10^ GeV this galactic component is not sufficient to describe 
the observed all-particle spectrum within the model. Hence, an additional component 
has to be introduced ad hoc in order to account for the measured flux values. Above 
3 • 10^ GeV the dashed line gives an upper limit for the proton fraction, assuming 
the ad-hoc component being protons only. The proton fraction decreases with energy 
from 35% at 10'* GeV to below 1% at 3 • 10^ GeV. Beyond this energy the upper limit 
increases again. However, taking the error band into account, at least between 10^ 
and 10* GeV protons contribute only little to the all-particle flux (< 5%). Thus, one 
of the premises to derive proton-air cross-sections in this energy region is weakened 
and, maybe, the values obtained from air shower measurements have to be corrected. 

In order to check the impact of inelastic cross-sections on the development of 
extensive air showers, the cross-sections are altered within a particular model. The 
present studies have been carried out using the simulation program CORSIKA version 
6.0190. Investigations by the KASCADE group (Antoni et al 1999, Milke et al 2001) 




Total inelastic cross- sections and X^ax in extensive air showers 



8 



40 



_l IIIIIU 1 IIIIIU 1 1 


iiiiij 1 iiiiiii| 1 iiiiuij 1 nil 


i| 1 1 null 1 1 Mill 


- 

— model 1 




- 








— model 3 








□ 


CDF : 




A 


E-710 - 




O 


E-811 - 


" 1 1 II 


Mill 1 IIIIIU 1 IIIIIU 1 MM 


1 1 1 null 1 1 Mill 


10^ 10^ 


5 6 7 8 9 
10 10 10 10 10 



Energy E^ [GeV] 



Figure 3. Inelastic proton-proton cross- 
section versus laboratory energy for three 
versions of the model QGSJET 01. The 
data points are for pp interactions from three 
experiments at the Tevatron, see text. 



indicate that the interaction model QGSJET 98 describes the measured data most 
consistently. Hence, the latest version of this code, i.e. QGSJET 01, has been 
chosen for the simulations. Within this program the paramctrizations for the inelastic 
hadronic cross-sections have been modified in order to change their energy dependence 
at high energies, the contribution of mini-jets has been reduced. The changes affect 
all inelastic cross-sections, like proton-proton, proton-air, and pion-air at the same 
time. In addition, a parameter had to be readjusted which essentially influences 
the multiplicity of produced particles in order to match pseudo-rapidity distributions 
of collider experiments. This implies a reduction of the average number of charged 
particles produced. 

In the following the original QGSJET code is referred to as model 1, while the 
two modifications with smaller cross-sections are labeled model 2 and 3, respectively. 
In addition, for model 3 the average transverse momentum has been reduced and the 
elasticity has been increased, these changes are combined in model 3a. 

The inelastic cross-sections for models 1 to 3 for proton-proton collisions are 
shown in figure 3 as function of laboratory energy. At 10^ GeV the cross-sections vary 
between a^J^' = 80 mb for the original QGSJET model and a™^' = 64 mb for model 3. 
At the Tevatron (i?o = 1-7 • 10® GeV) for the total elastic pp cross-section the values 
cr|^ 15.79 ± 0.87 mb (E-811, Avila et al 1999), cr|^ 16.6 ± 1.6 mb (E-710, Amos 
et al 1990), and af^ = 19.70 ± 0.85 mb (CDF, Abe et al 1994b) have been measured. 
Together with the total cross-sections mentioned above, values for the total inelastic 
cross-section are obtained as shown in the figure. Model 1 corresponds to tlic^ CDF 
value, model 2 to the results of E-710 and E-810, while model 3 marks approximately 
the 1(7 lower error bound for the last two experiments. 

Next, the model predictions are compared with air shower measurements. The 
proton-air cross-sections for models 1 to 3 have been plotted in figure 1 and are 
compared with experimental results. At 10^ GeV values between (J^^Iair — ™^ 
the original QGSJET 01 (model 1) and = 416 mb for model 3 arc obtained. 

At low energies up to 10^ GeV the results by Yodh et al (1983) and Mielke et al 
(1994) are compatible with the QGSJET cross-sections for all three cases considered. 
The cross-section obtained by the EAS-TOP collaboration is still lower, even when 
compared with model 3. On the other hand, the results from AGASA and Fly's 
Eye are about 1 to 2 cr above the calculated curve for model 3. Taking into account 
the low proton content of primary cosmic rays in this energy region according to the 
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( ). For comparison, cross-sections of interaction models implemented in 
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poly-gonato model, the experimental values are likely too large. The results of the 
EAS-TOP and AGASA experiments rely on the measurements of particle numbers at 
ground level. As has been pointed out by Alvarez-Muniz et al (2002a), intrinsic shower 
fluctuations strongly influence the cross-sections obtained from such measurements. 
The experimental uncertainties may be illustrated by two analyses of the Yakutsk data. 
The cross-sections obtained by Dyakonov et al (1990) are found to be compatible with 
model 3, while on the other hand, the more recent values of Knurenko et al (1999) are 
about 1(7 above model 3. Summarizing, no serious disagreement between experimental 
values and the cross-sections for model 3 can be stated. 

In figure 4 the cross-sections for models 1 to 3 are compared with results of other 
model calculations. In the left-hand panel the values for the models DPMJET 2.5, 
HDPM, SIBYLL 2.1, VENUS (Knapp et al 2003), and NEXUS 2 (Bossard et al 2001) 
as implemented in CORSIKA are shown. Compared with model 3 the cross-sections 
for DPMJET, NEXUS, and SIBYLL grow faster as function of energy. At energies 
below 10^ GeV the HDPM cross-sections are smaller than the values for model 3, 
while VENUS exhibits nearly the same behavior as model 3. 

The right-hand graph summarizes further theoretical considerations. Kopeliovich 
et al (1989) extrapolate results from collider experiments and calculate inelastic 
proton-air cross-sections. Their results for a QCD pomeron with an asymptotic 
intercept A = 0.097 are presented, these cross-sections are slightly larger than those for 
model 3. Also Frichter et al (1997) extrapolate results from accelerator measurements 
to high energies. They use relative low inelasticities and in turn obtain cross-sections 
larger as the previous model. But they admit that such low inelasticities are unable 
to account for the Fly's Eye data reported by Gaisser et al (1993) and Bird et al 
(1993). Similar values as for model 3 are obtained by Huang et al (2003) using the 
COSMOS simulation code. Gaisser et al (1993) investigate the falling slope of the 
observed X^^ax distribution as measured by the Fly's Eye experiment and conclude, 
that this decrement is consistent with a near linear lg(s) energy dependence of the 
inelastic proton-air cross-section. Their results are shown in the figure as well and are 
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Table 1. Inelastic proton-proton and proton-air cross-sections as function of 
laboratory energy for three versions of the high-energy hadronic interaction 
model QGSJET 01. Model 1 corresponds to the original QGSJET, in model 2 
and 3 the total inelastic cross-sections have been reduced, see text. 
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found to coincide with the values for model 3. The authors conclude further, when 
the extrapolations of proton-proton cross-sections from collider experiments to high 
energies by Block et al (1992) are converted to proton-air cross-sections (Gaisser et 
al 1987), one obtains values which are slightly above the lg(s) dependence mentioned 
above. Also these values are compatible with those of model 3. 

The inelastic cross-sections for both, proton-proton and proton-air interactions 
according to models 1 to 3 are listed in table 1 as function of energy for reference. 

3.2. Other model parameters 

The reduction of the contribution of mini-jets also influences the pseudo-rapidity 77. 
Therefore, another parameter mainly reducing the multiplicity had to be 
readjusted in order to match measurements at the Tevatron {^/s = 1.8 TeV ), shown 
in figure 5 together with the predictions of the models 1 to 3a. For central collisions 
(77 « 0) the values increase from model 3a to model 1. While for large jryj values the 
distributions for all models almost coincide. In conclusion, all four models yield very 
similar distributions, which describe well the CDF-measurements (Abe et al 1990), 
taking their error bars into account. 




-10 -8 -6 -4 -2 2 4 6 8 10 
pseudo rapidity r[ 



Figure 5. Pseudo-rapidity distribution of 
charged particles in p-p collisions at -^5 = 
1.8 TeV for four versions of QGSJET compared 
with Tevatron data (CDF, Abe et al 1990). 
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Figure 6. Number of charged particles produced in p-p collisions at =1.8 TeV 
for the QGS JET modifications together with measurements of E-735 (Alcxopoulos 
et al 1998) (left-liand side). Average rmmber of particles produced in tt-^^^N 
collisions as function of laboratory energy for four versions of the model QGSJET 
and values for three other models implemented in CORSIKA (Knapp et al 2003) 
(right-hand graph). 



The frequency distributions of the number of charged particles simulated for 
proton- antiproton collisions at ^/s = 1.8 TeV are plotted on the left-hand side of 
figure 6. The results obtained with the four models are very similar. A closer look 
reveals, the right-hand tail of the distributions is slightly steeper for models 3 and 
3a as compared with model 1. At the same time, the average multiplicity is reduced 
by up to 5% from 49.4 ± 0.1 (model 1) to 47.1 ± 0.1 for model 3 and 46.9 ± 0.1 for 
model 3a. The model predictions are compared with measurements at the Tevatron 
(E-735, Alexopoulos et al 1998). The data have been normalized to the calculated 
values for model 1 over a range of multiplicities starting at the peak of the distributions 
around Nch ~ 20. At medium multiplicities (w 100) model 1 slightly overestimates 
the production probability, while at high multiplicities all four models are below the 
measured values. 

Of great importance for the air shower development are pion-nitrogen collisions. 
The average multiplicity of charged particles (Nch) for these interactions is shown in 
figure 6 (right-hand side) as function of laboratory energy. For models 3 and 3a the 
average multiplicity is significantly lower as compared with model 1. The differences 
amount to more than a factor of three at the highest energies shown. Compared with 
other interaction models like DPMJET, NEXUS, or SIBYLL, relatively large values 
for (Nch) are obtained with the original QGSJET at high energies (Knapp et al 2003). 
On the other hand, the multiplicities obtained with models 3 and 3a coincide with the 
values for DPMJET 2.5. 

The average elasticity of the most energetic meson in pion-nitrogen interactions 
is plotted in figure 7 as function of laboratory energy for the four models. As can be 
seen in the figure, the reduction of the mini-jet contribution, as described above, also 
increases the elasticity at high energies. In model 3a the elasticity has been increased 
additionally by 10% to 15% relative to model 3 to check the influence of this parameter 
on the shower development. 

Reducing the inelastic cross-sections in QGSJET via a reduction of the mini-jet 
contribution also reduces the average transverse momentum {p±). Average values 
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Figure 7. Average elasticity of most energetic 

meson in tt- 

14n 

interactions versus laboratory 
energy for four versions of the model QGS JET. 



of p_L for charged pions and kaons in proton-antiproton collisions are shown in 
figure 8 (left-hand graph) as calculated with the different versions of QGSJET in 
the energy region of the Tcvatron. Measurements from experiment E-735 arc depicted 
for comparison. The average p±_ obtained with models 1 to 3 is up to 15% larger 
than the measurements for pions. Therefore, the average px predicted with model 3a 
has been adjusted to the experimental values shown. Unfortunately, the trend of the 
energy dependence of p±_ for kaons is less well described by all models discussed. But, 
since pions are clearly dominant for the air shower development, emphasis is given to 
a correct description of the transverse momenta for pions. 

The average transverse momentum of charged particles in pion-nitrogen collisions 
is shown in figure 8 (right-hand graph) as function of laboratory energy. At the highest 
energies the values for model 3a are about 15% below the average as calculated 
with model 1. 

In summary, it may be concluded that the four models exhibit distinct differences 

in certain quantities. Individual experimental results are described better by one or 
another model. But in a general view, all four versions of QGSJET are compatible with 




Figure 8. Left-hand panel: Average transverse momentum of pions (bottom) 
and kaons (top) in p-p collisions as function of cm. energy for four versions of the 
model QGSJET compared with Tcvatron data (E-735, Alcxopoulos et al 1993). 
Right-hand panel; Average transverse momentum of charged particles in 7r-^'*N 
interactions as function of laboratory energy. 
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results from collider and air shower experiments, or at least no serious disagreement 
between the models and the measurements can be stated. 

4. The average depth of the shower mELximum 

In this chapter, the consequences of changing the inelastic cross-sections and the 
other modified parameters on the longitudinal development of air showers shall be 
investigated. For this purpose, in the energy range from 10^ GeV up to 3 • 10^° GeV, 
500 showers induced by primary protons and 100 by iron nuclei have been simulated 
for each energy in steps of half a decade using the CORSIKA program. To describe 
the high-energy hadronic interactions the models 1 to 3a have been used together with 
the GHEISHA code for energies below 80 GeV. The discussion of the results is divided 
in two sections, on average Xj^ax values and on X^^x distributions. 

Average X^^ax -values 

The average depth of the shower maximum as obtained by the calculations is shown in 
figure 9 as function of energy for primary protons and iron nuclei for the four models 
considered. The original QGSJET yields the smallest Xmax values, with model 3a the 
showers penetrate deepest into the atmosphere. The values change at 10^ GeV from 
Xmax = 668 ± 3 g/cm2 to X 

max — 702 ± 3 g/ cm^ for proton induced showers and 




Figure 9. Average depth of the shower maximum versus primary energy. Results 
are shown from the experiments BLANCA (Fowler et al 2001), CACTI (Paling 
ei al 1997), DICE (Swordy et al 2000), Fly's Eye (Bird et al 1994), Haverah 
Park (Watson 2000), HEGRA (Arqueros et al 2000), HiRes/MIA (Abu-Zayyad 
et al 2000), Mt. Lian Wang (Cha et al 2001), SPASE/VULCAN (Diekinson et al 
1999), and Yakutsk (Dyakonov et al 1993, Kimrenko et al 2001). The data are 
compared with results from CORSIKA simulations for primary protons and iron 

nuclei using different versions of QGSJET: model 1 (original QGSJET) ( ), 

model 2 ( ), model 3 ( — • — ), and model 3a ( ). 
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Table 2. Average depth of the shower maximum for primary protons and iron 
nuclei calculated with CORSIKA, using different versions of the interaction 
model QGSJET, sec text. 



Average depth of the shower maximum Xmax [g/cm^] 
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from X„iax = 567 ± 2 g/cm^ to X„iax = 577 ± 2 g/cm^ for iron primaries. X^ax is 
listed in tabic 2 as function of energy for the four versions of the QGSJET model. 

The simulated values are compared with experimental results from various 
experiments in the figure. Two different techniques are used to measure Xmax, namely 
the imaging technique, using telescopes to obtain a direct image of the shower and 
the non-imaging method, in which the height of the shower maximum is derived 
from the lateral distribution of the Ccrenkov light (or electrons in case of Haverah 
Park) measured at ground level. The DICE, Fly's Eye, and HiRes experiments use 
the imaging method, while all other experiments belong to the second group. The 
data show systematic differences of « 30 g/cm^ at 1 PeV increasing to ^ 65 g/cm^ 
close to 10 PeV. A more detailed discussion of the measurements is given elsewhere 
(Horandel 2003). Some of the experimental uncertainties may be caused by changing 
atmospheric conditions. The experiments measure a geometrical height which has to 
be converted into an atmospheric depth. Measuring longitudinal atmospheric profiles 
during different seasons, Keilhauer et al (2003) found that the atmospheric overburden 
for a fixed geometrical height (e.g. 8 km a.s.l.) varies by at least 25 g/cm^. 

The ambiguities between the different experiments are of the same order as the 
systematic differences between individual models as discussed in chapter 2, where 
ambiguities of about 45 g/cm^ for proton and about 25 g/cvo? for iron induced 
showers have shown up. It can be noticed that up to 4 PeV the measured Xmax values 
increase faster as function of energy than the simulated values for the original QGSJET 
(model 1). With the modifications in model 3a the sinndated showers penetrate deeper 
into the atmosphere and the increase is almost parallel to the data. 

An estimate for the dependence of JS.Xmax/ X^ax on the multiplicity and the 
inelasticity has been given in equation (1). The present investigations confirm the 
dependence of Xmax on the inelasticity. The correlations between changes in the 
inelastic proton-proton as well as proton-air cross-sections and Xmax can be estimated 



Total inelastic cross- sections and X^ax in extensive air showers 15 



by 

/\Y 9 Arr""'' \Y Arr*"^' . 

^^max ^ ^ ^^-^ pp Qxid ^ p—air /^x 

-^max ' ^pp -^max ' ^p—air 

Since the changes of the parameters in QGSJET affect aU cross-sections 
simuhaneously, the dependencies on one particular cross-section can not be 
disentangled. 

In the superposition model of air showers, the development of a cascack; induced 
by a primary particle with energy E and nuclear mass number A is described by A 
proton induced subshowers of energy E/A. Hence, the shower maximum of a proton 
induced shower X^^^ should be equal to that for an iron primary X^^^ at higher 
energy 

XP^,^{E) ^ X^l^{E ■ A^e) . (4) 

The quantity ^{E) defined as 

1^{E) EE X^l,{E ■ Ap,) - X^l,{E) = Xl,,{E) - Xil,{E) (5) 

specifies the increase of the Xmax curve in the energy interval [E, Ape ■ E), it relates 
to the frequently used elongation rate V = dX^ax/dlgE as A{E) = 'D{E)\g{Ape). 
The right-hand side of equation (5) follows immediately from equation (4). This 
indicates, that the slopes of the Xmax curves are correlated with the difference 
between the proton and iron curves at a given energy. Such a behavior is visible 
in figure 9. For model 3a the curves for protons and iron nuclei are almost straight 
lines in the logarithmic plot. In the energy range shown, A{E) changes from about 
150 g/cm^ to approximately 120 g/cm^ only. On the other hand, the logarithmic 
slope of the two curves for model 1 decreases significantly with rising energy and the 
difference between protons and iron diminishes from A(3 • 10^ GeV) ~ 150 g/cm^ to 
A(3 • 10^° GeV) « 80 g/cm^. 

In the previous chapter cross-sections for several models have been presented. It is 
interesting to compare the depths of the shower maxima calculated with these codes 
to the depths obtained with the QGSJET modifications. Xmax values for primary 
protons and iron nuclei for the interaction models DPMJET 2.5, NEXUS 2, and 
SIBYLL 2.1 have been calculated by Knapp et al (2003) using the CORSIKA program. 
The results are compared in figure 10 to values for models 1 to 3a. SIBYLL predicts 
relatively early developing cascades with values between models 1 and 2. On the other 
hand, DPMJET generated showers penetrate deeper into the atmosphere similar to 
model 3a. The validity of the NEXUS code is guaranteed by its authors up to 10*^ GeV 
and, correspondingly, Xmax up to this energy is plotted. The depths calculated 
with this code are very similar to the SIBYLL predictions. Depths simulated with 
the MOCCA code using the internal interaction model (Prykc 2001) are significantly 
larger as compared with the other codes. With the original QGSJET 01 (model 1) 
the shortest development of the cascades is obtained. 

In figure 4 it has been shown, that the inelastic cross-sections for SIBYLL rise 
faster as compared with DPMJET. This is compatible with deeper penetrating showers 
for DPMJET, as can be inferred from figure 10. How also the other parameters, like 
the multiplicity, effect the longitudinal development can be seen as an example with 
DPMJET. Its showers penetrate deeper into the atmosphere than those generated 
with model 1, despite the almost identical cross-sections (see figure 4). However, the 
multiplicity of charged particles grows faster with energy for QGSJET as compared 
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Figure 10. Average depth of shower 
maximum versus primary energy. Predictions 
of simulations for proton and iron induced 
showers using CORSIKA with the interaction 
models DPMJET, NEXUS, and SIBYLL 
(Knapp et al 2003) as well as the MOCCA 
program (Pryke 2001) are compared with 
results of the modifications of QGSJET 01 
(models 1 to 3a). 

with DPMJET, e.g. at 10^ GeV the multiphcities differ by about a factor of 3 (see 
figure 6). 

4-2. X max -distributions 

In addition to the average depths, probabihty distribiitions of X^ax values have 
been pubhshed and can be investigated. In the energy range from 3 • 10^ to beyond 
10^ GeV experimental Xma^-distributions from the Fly's Eye (Gaisser et al 1993) and 
HiRes (Abu-Zayyad et al 2001) experiments are presented in figure 11. Early Xmax 
distributions of the Yakutsk experiment (Dyakonov et al 1990) are not compatible 
with more recent average Xmax values (Dyakonov et al 1993), shown in figure 9, and 
are therefore not taken into account. 

The results of the measurements are compared with calculated Xmax distributions 
using model 3a for primary protons, helium, oxygen, iron, and ultra-heavy nuclei. 
The CORSIKA code allows primary nuclei with mass numbers A < 56. The Xmax 
distributions for ultra-heavy nuclei have been estimated assuming Xmax oc In^, based 
on the simulated values for protons and iron nuclei. For each energy interval the mass 
composition is taken according to the poly-gonato model, i.e. the abundances are taken 
as predicted and listed in table 3. Distributions for five mass groups, represented by 
the nuclei mentioned above, have been calculated. The resulting X^ax-distributions 
are shown in figure 11 as solid histograms. For the three energy intervals the left- 
hand side of the Fly's Eye distributions is described reasonably well by the galactic 
component. It should be stressed that the agreement has not been achieved by a fit 
to the distributions. Instead, the mass composition as obtained with the poly-gonato 
model has been chosen. 

When comparing the same Fly's Eye data to results obtained with the original 
QGSJET code, an offset AXmax = 30 g/cm^ had to be added to the simulated Xmax 
values in order to reconcile the measured with the simulated distributions (Horandel 
2003). A shift in the opposite direction has been performed by Gaisser et al (1993) in 
order to achieve agreement between calculations with the KNP model and the Fly's 
Eye data, the simulated values have been shifted by AXmax = —25 g/cm^. As can 
bee seen in figure 11, with model 3a no artificial offset is necessary to describe the 
Fly's Eye distributions at the side of low Xmax- This indicates experimental evidence 
for deeper penetrating showers as obtained with model 3a. 

The galactic component accounts only for a fraction of the observed cosmic rays 
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Figure 11. Distribution of the depth of the 

shower maximum Xmax measured by the Fly's 
Eye (Gaisscr et al 1993) and HiRes (Abu- 
Zayyad et al 2001) experiments. The measured 
values are compared witli simulated results 
using model 3a. Galactic component according 

to the poJy-gonato model ( ) and galactic 

plus ad-hoc component ( ), see text. 



in the energy range considered in figure 11. As mentioned in chapter 3, close to 
10® GeV the heaviest elements of the galactic component reach their cut-off energies 
and a new component starts to dominate the all-particle spectrum. For the additional 
component only the flux but not the mass composition can be predicted. Assuming 
only protons for this component cannot explain the experimental X^ax distributions. 
Hence, helium and the CNO group have been added. A fit to the Fly's Eye data 
using the simulated distributions for protons. hcliTim, and oxygen yields the dashed 
histograms as the sum of the galactic and the ad-hoc component. The relative 
abundances obtained are listed in table 3 in brackets. As can be seen in the figure, 
the dashed histograms reproduce well the Fly's Eye distributions, indicating a mixed 
mass composition for the ad-hoc component. The contribution of the CNO group 
seems to be essential to describe the peak region of the experimental distributions. 
Protons account for the long, asymmetric tail towards large values of Xmax- The steep 
gradients on the left-hand side originate in the relative narrow distributions for the 
ultra heavy elements of the galactic component with their small intrinsic fluctuations 
in the shower development. 

Using the HiRes data for the fit would result in a slightly lighter mass composition. 
But due to the relatively large errors for the energy range 5 — 10 • 10® GeV a fit seems 
not to be meaningful. 

The fraction of protons obtained for the ad-hoc component is plotted in figure 2 
as function of energy (filled points). Between 10^ and 10^ GcV the primary cosmic-ray 
fiux is not dominated by protons. Between 6-10^ GeV and 5-10® GeV the fraction 
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Table 3. Relative abundances [%] for groups of nuclei with charge number Z 
for different primary energies Eq for the galactic component according to the 
poly-gonato model. Values in brackets are for the ad-hoc component in order to 
explain X-max measurements of the Fly's Eye experiment. 
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amounts to less than 10%, between lO'' GeV and 2 • 10^ GeV its value is even below 
5%. One has to bear in mind this indication when proton-air cross-sections are derived 
from air shower measurements. 

The low inelastic cross-sections and higher values for the elasticity assumed in 
model 3a have implications on the mass composition derived from Xmax measurements 
at energies above 10'' GeV as pointed out above. With the low cross-sections also 
intermediate and heavy elements are important in this energy region. In the energy 
region between 10'' and 10* GeV a relatively heavy composition has been found also 
by Alvarez-Mufiiz et al (2002b), viz. consisting of 85% Fe, 10% CNO, 4% He, and 
1% protons. Similar values are calculated with the poly-gonato model at 4 • 10'' GeV: 
86% heavy and ultra-heavy (9 < Z < 92), 9% CNO (6 < Z < 8), 4% He (2 < Z < 5), 
and < 1% protons. A key issue of the present investigations is that cosmic rays above 
10* GeV, i.e. the ad-hoc component, contain a significant contribution of particles 
heavier than protons. A mass composition heavier than protons only in this region is 
also obtained by Erlykin and Wolfendale (2002). 



5. Mean logarithmic mass 

Many scientists characterize the mass composition of high-energy cosmic rays by the 
mean logarithmic mass. It is defined as 

{\YiA) = ^r,\nA, , (6) 

i 

with the relative fraction of nuclei with mass Ai. Knowing the average depth of the 
shower maximum for protons X'p^^^^ and iron nuclei X^^^ from simulations, the mean 
logarithmic mass can be derived in the superposition model from the measured values 
XTaT using 

vmeas vp 

{lnA)= — _ — -In^^e . (7) 

The corresponding (In A) values for the variations of QGSJET 01, obtained from 
the data presented in figure 9, are plotted versus the primary energy in figure 12 for 
models 1 and 2 as well as for models 3 and 3a in figure 13. The average (In A) increases 
as the cross-sections decrease from model 1 to model 3. For the original QGSJET the 
results of many experiments exhibit a (strong) decrease of (In A) up to about 4-10^ GeV 
and an increase above this energy. The energy of this dip in the (In ^4) values coincides 
with the energy of the knee in the all-particle energy spectrum. The dip becomes less 
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Figure 12. Mean logarithmic mass (in A) of primary cosmic rays versus their 
energy calculated from experimental Xmax values from air shower observations 
(see figure 9) with the results of CORSIKA/QGSJET simulations for two different 
sets of inelastic hadronic cross-sections, model 1 (top) and model 2 (bottom). 
Results of direct measurements from JACEE (Shibata 1999) and RUNJOB 
(Apanasenko et al 2001) are shown as well. The lines represent calculations 

according to the poly-gonato model: galactic component only ( ), plus ad-hoc 

component of solely protons ( ). 



striking with lower inelastic cross-sections and higher values for the elasticity. For 
model 3a only a modest dip can be recognized. At 4 PeV the average values increase 
from (In A) = 1.2 for model 1 to (In A) = 1.6 for model 3a. Around 10** GeV the 
average logarithmic mass compared with model 1 is about A{\nA) « 0.5 larger for 
model 3 and A(ln^) ^ 0.7 larger for model 3a. These examples illustrate how strong 
the interpretation of air shower measurements depends on model parameters like the 
inelastic cross-sections or elasticities used. At Tevatron energies the cross-sections 
vary within the error range given by the experiments and at 10^ GeV the proton-air 



Total inelastic cross- sections and X„iax in extensive air showers 



20 




Figure 13. Mean logarithmic mass (in A) of primary cosmic rays versus their 
energy calculated from experimental Xmax values from air shower observations 
(see figure 9) with the results of CORSIKA/QGSJET simulations using models 3 
(top) and 3a (bottom). See also caption of figure 12. 



cross-sections of models 1 and 3 differ only by about 10%, but the general trend of 
the emerging (In^) distributions proves to be significantly different. 

At this point the circle closes. If we assume in the energy region from 10^ to 
10® GeV a small proton fraction only, the cross-sections have to be corrected and 
lowered to the region of values for model 3. In turn a heavier composition is obtained 
for model 3 as has been demonstrated in figure 13. Thus, at least qualitatively the 
arguments are consistent. 

Results from the balloon experiments JACEE (Shibata 1999) and RUNJOB 
(Apanasenko et al 2001) are presented in figures 12 to 14 (filled points) for comparison. 
No hint for a decreasing mean logarithmic mass is indicated by these measurements. 
The solid lines shown in the figures are predictions according to the poly-gonato model 
for the galactic component and the dashed lines include an ad-hoc component of 




Figure 14. Mean logarithmic mass (In A) of primary cosmic rays versus their 
energy calculated from experimental Xmax values from air shower observations 
(see figure 9) with the results of CORSIKA simulations for three different hadronic 
interaction models as indicated, see figure 10. See also caption of figure 12. 
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protons only. Both, the results of the balloon experiments shown as well as the values 
calculated with the poly-gonato model seem to support model 3a. 

In addition to the modifications of QGSJET, it would be interesting to look 
at the mean mass which is deduced from the measurements applying other models. 
Xmax values for the models DPMJET, NEXUS, and SIBYLL have been presented in 
figure 10. Taking these predictions the experimental results lead to mean logarithmic 
masses shown in figure 14. A t 4 PeV the average mass values range from (In A) = 1.7 
for SIBYLL and NEXUS to (In^) = 2.0 for DPMJET. For NEXUS the authors of 
the model guarantee validity up to 10* GeV, as mentioned above, for this reason 
(In A) is shown up to this energy only. As can be inferred from figure 10, very similar 
depths arc obtained for the models NEXUS and SIBYLL, both are hardly discernable 
from model 2. As a consequence, the mean logarithmic masses derived are very much 
alike. On the other hand, the depths predicted by DPMJET are very similar to 
the results of model 3a and for both models a heavier mass composition is obtained. 
For the program MOCCA with its internal interaction model the showers penetrate 
extremely deep into the atmosphere and a relative heavy mass composition is obtained, 
not compatible with direct measurements (Horandel 2003). 

At knee energies the average experimental (In^) values vary from 1.2 when the 
data are interpreted with the original QGSJET (model 1) to 2.0 for a DPMJET 
interpretation. One has to admit that the model ambiguities result in an uncertainty 
A(ln A) « 0.8. At this energy the scatter in the measured average depth of the shower 
maximum, as presented in figure 9, yields a rms value of the individual experimental 
results in the order of rms^jn^j w 0.6. This value is compatible with rms^in^^ w 0.6 
obtained for experiments measuring particle distributions at ground level, as discussed 
elsewhere (Horandel 2003). 

Also other authors have studied the effects of different interaction models on 
(In A). The systematic influence of the models HDPM, QGSJET, SIBYLL, and 
VENUS on the results of the BLANCA experiment is discussed by Fowler et al 
(2001). At knee energies a maximum offset A{lnA) » 0.8 between QGSJET and 
SIBYLL as well as HDPM is obtained. The investigations of Wibig (2001) yield a 
maximum difference A(ln A) « 0.9 between VENUS and SIBYLL. Comparing results 
from QGSJET 98 and SIBYLL 1.6 an uncertainty of A(ln A) » 0.3 is foimd by Swordy 
et al (2002) at knee energies. Investigations of several hadronic observables by the 
KASCADE group yield model ambiguities of A(ln A) w 0.4 around the Jaiee (Horandel 
et al 1998). These estimates from the literature seem to be well compatible with the 
ambiguities determined above. 

As already mentioned in the introduction, the mean logarithmic mass as obtained 
with experiments investigating particle distributions at ground level is compatible with 
the results of the poly-gonato model. But these results disagree with {In A) values 
obtained from measurements of the longitudinal shower development interpreted with 
QGSJET 01. This incompatibility can be seen in figure 12, as well as for the models 
NEXUS 2 and SIBYLL 2.1 presented in figure 14. For several experiments the mean 
logarithmic mass decreases as function of energy up to knee energies, a tendency not 
supported by the extrapolation of the direct measurements. This effect is strongest 
for QGSJET 01. 

Using lower cross-sections and larger values for the elasticity in the model 

QGSJET the {In A) values obtained are comparable with the results of experiments 
investigating particle distributions. In other words, the disagreement between the 
two groups of experiments can be reduced, if model 3a is taken. With the altered 
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inelastic cross-sections and the larger elasticity consistency can be achieved between 
the predictions of the poly-gonato model and the mass composition derived from 
observed X„iax values. Similar values of (In A) as for model 3a are also obtained using 
DPMJET 2.5 to interpret the data. The main conclusion of the present investigation is 
that relatively deep penetrating showers with X^ax values similar to the ones obtained 
with model 3a or DPMJET 2.5 seem to yield most consistent (in A) values. On the 
other hand, investigations of secondary particles produced in air showers by Milke et al 
(2001), based on CORSIKA simulations with the low energy model GHEISHA, reveal 
that for a given number of muons DPMJET 2.5 transports more hadronic energy to the 
observation level as compared with the measurements of the KASCADE experiment. 
Hence, presently model 3a is the preferred model to describe air shower measurements 
most consistently. 

6. Number of electrons and muons 

The modified cross-sections and elasticity values influence not only the average depth 
of the shower maximum but also other air shower observables. The implications on 
shower sizes at ground level might be strong and, therefore, are investigated in the 
following. Many air shower arrays use the correlation between the number of electrons 
Ne and the number of muons Nf^ to determine the cosmic-ray mass composition. 
For showers with primary energies between 10^ and 10^*^ GeV the average number 
of muons with energies above 100 MeV is plotted versus the average number of 
electrons with energies above 0.25 MeV in figure 15 (left-hand side). Correlations 
from CORSIKA/QGSJET simulations for proton and iron induced showers are shown, 
using the original QGSJET 01 and model 3a. The first impression is, no significant 
differences between the two models can be found and, consequently, no significant 
changes in the mass composition derived are expected. 

A closer look, when investigating the relative changes for the numbers obtained 
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Figure 15. Left-hand panel: Average number of muons (E^ > 100 MeV) versus 
number of electrons (Ee > 0.25 MeV) for model 1 and 3a. The range pictured 
corresponds to primary energies from 10^ to 3 ■ lO^" GeV. Right-hand panel: 
Number of muons and electrons as predicted by model 3a relative to the numbers 
calculated with model 1 as function of primary energy. The range for the ratio 
Af^"/Af^ for the extrema protons and iron nuclei is illustrated by solid lines and 
the range for the ratio N^'^ /N^ is represented as shaded area. 
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reveals small deviations. The number of electrons and muons obtained with model 3a 

normalized to the results of the original QGSJET model are shown on the right- 
hand side of the figure. At energies of the knee (4 PeV) the number of electrons 
increases by about 5% when using model 3a for both primaries, protons and iron nuclei. 
The increase in the number of muons amounts to about 15%. This outcome may 
qualitatively explain the need to accommodate the energy scales of many air shower 
experiments with respect to direct measurements as outlined in the introduction. It 
has been found that at knee energies air shower experiments overestimate the energy 
on average by about 3%. This can be understood, if the interaction models predict 
too small electron and muon numbers, since these quantities are mostly used by 
experiments to estimate the primary energy. In fact, the change in electron and 
muon numbers from models 1 to 3a corresponds to an energy uncertainty of about 5% 
and 15%, respectively. These numbers depend on the threshold for electrons as well 
as muons and differ from experiment to experiment. A quantitative investigation for 
each experiment would require detailed detector simulations. However, the numbers 
indicate that the observed energy overestimation of air shower experiments relative 
to direct measurements may be explained by too high inelastic cross-sections and too 
low values for the elasticity. 

At higher energies around 10^ GeV the number of both, electrons and muons, 
generated by model 3a is increased by about 15% to 20% as compared with the original 
QGSJET (model 1). The changes in A^e and correspond to energy shifts in the 
order of 10% to 20%. In this energy region the air shower experiments overestimate 
the primary energy on average by about 10% according to the above mentioned 
investigations. This value is comparable with the systematic offset between models 1 
and 3a. Again, the exact numbers depend on individual detection thresholds for each 
experiment. These findings are an independent hint for a slower logarithmic rise of 
the inelastic cross-sections as well as slightly increased elasticities and, consequently, 
deeper penetrating cascades as presumed so far. 

It may be summarized that differences between the models shift the data points 
parallel to the lines on the left-hand side of figure 15 (i.e. in "energy direction") and 
not perpendicular (i.e. in "mass direction"), hence only a change in energy and no 
significant change in the mass composition derived from the data are expected. More 
detailed investigations are necessary to study the implications of the modifications on 
observables, performing detector simulations for individual experiments, but this is 
beyond the scope of the present article. In addition, to interpret experimental data 
not only high-energy interaction models as discussed presently are important, but also 
models to describe interactions below 100 GeV, as for example GHEISHA, FLUKA, 
or UrQMD, all available in CORSIKA. For a final conclusion, the complex interplay 
between low and high-energy models as well as their influence on observables measured 
in ground arrays has to be studied. 

7. Summary and conclusion 

The impact of reduced inelastic cross-sections and increased values for the elasticity 
on the development of extensive air showers has been studied with the simulation 
program CORSIKA. Within the high-energy hadronic interaction model QGSJET 
the logarithmic increase of the inelastic cross-sections as function of energy has 
been lowered and the elasticity has been increased. The average depth of the 
shower maximum was calculated for primary protons and iron nuclei and compared 
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with experimental values. The mean logarithmic mass (In A) was derived from the 
experimental data for three assumptions for the increase of the inelastic cross-sections 
and two hypotheses for the elasticity. 

The disagreement in the observed (In A) values between experiments measuring 
particle distributions at ground level and experiments measuring the average depth 
of the shower maximum could be reduced. The data exhibit still a scattering in the 
order of A(lnyl) « 0.5 but the general trend of the increase as function of energy is 
now similar for both classes. A reduced and only modest increase of the cross-sections 
as function of energy had to be applied. Best agreement is obtained for model 3a with 
an increase of the total inelastic proton-proton cross-section to 72 mb at 10^ GeV and 
an additional increase of the elasticity of about 10% to 15%. In turn, the increase of 
(In A) with energy according to the poly-gonato model becomes compatible with the 
two classes of experiments. As a result, model 3a allows a consistent description of the 
extrapolations of individual element spectra as obtained by direct measurements and 
the all-particle energy spectrum as well as the cosmic-ray mass composition obtained 
by most air shower experiments. 
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